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ABSTRACT 

We present the results of timing analysis of the low-frequency Quasi-Periodic Oscillation (QPO) 
in the Rossi X-Ray Timing Explorer data of the black hole binary XTE J1550-564 during its 1998 
outburst. The QPO frequency is observed to vary on timescales between '^100 s and days, correlated 
with the count rate contribution from the optically thick accretion disk: we studied this correlation 
and discuss its influence on the QPO width. In all observations, the quality factors (I'o/FWHM) 
of the fundamental and second harmonic peaks were observed to be consistent, suggesting that the 
quasi-periodic nature of the oscillation is due to frequency modulation. In addition to the QPO and 
its harmonic peaks, a new l.Biy component was detected in the power spectra. This component is 
broad, with a quality factor of ~0.6. From this, we argue what the peak observed at half the QPO 
frequency, usually referred to as "sub-harmonic" could be the fundamental frequency, leading to the 
sequence 1:2:3:4. We also studied the energy dependence of the timing features and conclude that 
the two continuum components observed in the power spectrum, although both more intense at high 
energies, show a different dependence on energy. At low energies, the lowest-frequency component 
dominates, while at high energies the higher- frequency one has a higher fractional rms. An interplay 
between these two components was also observed as a function of their characteristic frequency. In 
this source, the transition between low/hard state and hard-intermediate state appears to be a smooth 
process. 

Subject headings: accretion: accretion discs — black hole: physics — stars: individual (XTE J1550- 
564) — X-ray: binaries 



1. INTRODUCTION 

Since the launch of Rossi X-ray Timing Explorer 
(RXTE), an extraordinary progress has been achieved 
in the knowledge of the variability properties of black- 
hole candidates (BHCs) in X-ray binaries. Different 
types of Quasi-Periodic Oscillations (QPO) have been 
observed in these systems. While only a few bina- 
ries show high-frequency QPOs (HFQPOs, 50-450 Hz , 
see iRemillard fc McClintoclj [20061 : iBelloni et all l2006l ). 
low-frequ ency QPOs (LFQPOs, mHz to ^10 Hz , see 
[RemiUard fc McClintockl [20061: ICaseTla et all 120051) are 
detected in virtually all observed BHCs. As both low 
and high-frequency QPOs are thought to arise in the ac- 
cretion flow close to the black hole, the study of their 
properties and behavior can provide important clues on 
the physics of accretion onto BHCs. 

In the case of LFQPOs, several distinct types show- 
ing different properties have been observed. Three main 
types, dubbed types A, B, and C respectively, s tand 
out in the pr esent scenario. iWiinands et~an (|1999[ ) and 
iHoman et al.l (2001) first report ed type A and type B 
LFQPO in XTE J1550-564, while lRemillard et al.l (|2002[ ) 
dubbed the ubiquitous LFQPO that appears together 
with band-limited noise type-C. The detailed properties 
of the different types of LFQPOs were investigated by 
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iCasella et al.l ()2005[ ). However, the physical difference 
among types remains unknown. 

The X-ray transient X TE J1550-564 was discovered 
on 1998 September 7th (iSmith et al.l I1998D with the 
RXTE All Sky Monitor (IWood et al.lll9990 . The discov- 
ery prompted a follow-up series of almost daily pointed 
RXT E/PCA (Proportional Counter Arrav: [Jahoda et al.l 
Il996f ) observations, which revealed a hard power-law 
dominated spectrum. Two weeks later, it reached a peak 
intensity of 6.8 Crab at 2-10 keV. The marked softening 
of the spectrum during this period indicates a transi- 
tion from low/hard state (LS) to hard-intermediate state 
(HIMS), reaching the soft-intermediat e state ( SIMS) in 
occasion of the 6 .8 Crab peak (see IHoman fc Bellonil 
[20051: lBellonll2009l ). After the bright peak, XTE J1550- 
564 remained in the HIMS for more than three weeks. 
Strong LFQPOs were observed, with frequency changing 
quite dramatically during the first ~ 40 d of the out- 
burst but then stabilizing from day 40-52. The type 
of QPO also changes at day 40 (from C/C to B: see 
IRemillard et al.ll2002D . Additional outbursts of the sys- 



tems followed in 2000, 2002 and 2003. 
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Subsequent optical observations showed that the dy- 
namical mass of the compact object is 10.5 ± 1.0 Mq, 
indicating a black-hole nature. Its binary companion 
was found to be a low-mass star, and the dis tance to the 
sourc e w as estimated t o be a bout 5.3 kpc (jOrosz et al.l 
I2002f) . iCorbel et all (|2002D discovered a large-scale, 
relativistically moving and decelerating jet emitting in 
radio and X-rays. 
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In this paper, we concentrate on the power density 
spectra found at the beginning of the 1998 outburst of 
XTE J1550-564 and analyze in detail the type-C QPO 
and noise components, focussing on their relative prop- 
erties. Particular attention is given to the "harmonic" 
peaks, which reveal information about the nature of the 
observed signal. 

2. OBSERVATIONS AND DATA ANALYSIS 

IRemillard et al.l (|2002[ ) analyzed the X-ray power den- 
sity spectra (PDS) for ah the 209 RXTE observations of 
XTE J1550-564 during its major outburst of 1998-1999, 
and detected all three types of LFQPOs. Most of them 
are of type C and occu rred in the first half of the outburst 
(jSobczak et al.|[20Q0bf ). We reanalyzed the data and fo- 
cussed on the 47 observations with type-C QPO between 
MJD 51,065 and MJD 51,101. We also analyzed ah nine 
observations with type-B LFQPO for comparison. 

Custom timing-analysis software under IDL and MAT- 
LAB was used. For each observation (In a few cases, 
an inspection of the PDS showed significant variations 
in the QPO frequency between different RXTE orbits, 
which were therefore split.), we produced a set of PDS 
from 128-s long stretches from the PCA channel band 
0-35 (corresponding to 2-13 keV). The time resolution 
was 1/256 s, corresponding to a Nyquist frequency of 
128 Hz. These spectra are then averaged together and 
logarithmically rebinned. The effect of dead time on 
Poisson noise was not subtracted directly, but fitted in 
the power spectra with an additional additive constant. 
Broad-band and peaked featur es wer e modeled by a com- 
bination of Lorentzians (|Nowakll2000l : iBelloni et aIll2002D 
using XSPEC vl l.3. The power sp ectra were normal- 
ized according to iLeahv et al.l (|1983l ): however, the fit- 
ting results were convert ed to squared fractional rms 
(|BeUoni fc Hasinged[T990l ). 

For each observation, a spectrogram was also accumu- 
lated, in order to detect time variations of the QPO pa- 
rameters within a single observation. This consists in a 
time sequence of the single 128-s power spectra, i.e. a 
time-frequency image. We fitted each PDS in the spec- 
trogram (limited to a narrow frequency range centered 
on the main QPO peak) with a model consisting of a 
Lorentzian peak and a power law for the local contin- 
uum. In this way, we derived the frequency shift with 
time on a 128-s time scale. 

For each observation, we also produced PDS in differ- 
ent channel ranges: eight sub-bands of the 0-35 range 
plus channels 36-89 (corresponding to 13-33 keV). 

2.1. PDS model 

In order to characterize the LFQPO behavior, a con- 
sistent model is required to describe the full PDS for 
the observations showing a type-C QPO. During the 
initial rise of the outburst, a state transition from LS 
into HIMS takes place: the spectrum softens considcr- 
ably and the Q PO frequencies increase well abov e 1 Hz 
(iCui et al.l [1999 : Sobczak et_ al IQQOb; Rcmilla rd et al.l 
|2002| ). Therefore, changes are expected in the overall 
shape of the PDS. We find that the PDS of first 5 ob- 
servations appear different (see Fig. [T|). For them, we 
adopted a model consisting of a fiat-top noise component 
L/t, one or two band-limited noise components Lblni/2^ 



a QPO peak Lp with its second harmonic QPO Lh (see 
Fig. [l] panel a). For all other observations with type- 
C QPOs, we used a model consisting of a fiat-top noise 
component L/t, a peaked noise component Lp„, a QPO 
Lp, with a sub-harmonic and a second harmonic Lh 
(see Fig. [H panel b). For some PDS a third harmonic 
appears. These two models fit the data reasonably well, 
with best-fit reduced values less than 2 (for ~ 265 
dof), with a typical value of 1.5. Hereafter, we con- 
centrate on the second part of observations, which we 
identify as a HIMS. As presented below, the results of 
our analysis suggest that the fundamental frequency of 
the oscillation is Lg rather than the conventional Lp- 
However, for clarity, we will refer to Lp as the funda- 
mental throughout the paper. For all noise and QPO 
components, we consider their characteristic frequency 

Vmax = \/ + (^/2)^, where and A are the centroid 
frequency and the FWHM of t he Lorentzian peak respec- 
tively (see IBelloni et al.l 120021 for a discussion). This is 
the frequency at which the vPu power spectrum peaks. 
In the case of broad components, this is a more ratio- 
nal choice than the centroid frequency, since the use of 
Lor entzian mode l s does not have a p hysical motivation 
(see lNowak|[2000t IBelloni et al.l [2003) . With this defini- 
tion, it is possible to compare homogeneously narrow and 
broad features, which have been to evolve from one to the 
other (see e.g. iDi Salvo et all[200l[ ). and to discover ma- 
jor correlations between characteristic fre quencies both 
in neutron-star and blac k-hole systems (see lBelloni et al.) 
[200l [ van der KlisI 120061 ). The frequencies of component 
X, Lx, will be indicated as vx- 

For the PDS with a type-B QPO, we used a model 
consisting of a power-law for red noise, two Lorentzians 
for the sub-harmonic QPO and the second harmonic 
QPO, plus a Gaussian for the fundam ental QPO (see 
iNespoh et alll2003t [Casella et al.ll2004D . The best-fit re- 
duced values were similar to those obtained for type-C 
power spectra. 

3. RESULTS 

3.1. Correlations between frequencies 

Following IBelloni et al.l (|2002D . we calculated the char- 
acteristic frequency Vmax of each component in the PDS 
of the HIMS observations. In Fig. [2l we plot all the fre- 
quencies versus that of the QPO [Lp). It is clear that 
the frequency ratio of the three QPO peaks Vg, up and 
Vh is 1/2:1:2, confirming their harmonic relation. The 
frequency of the flat-top noise component Vft appears to 
to be a factor of ~ 5 below vp, consistent with the c orre- 
lation repo rted by Wii nands fc van der KlisI ([19991) and 
IBelloni eta l. (2002). Interestingly, the characteristic fre- 
quency of the peaked noise component fp„ is ~ 1.5 times 
up^ although some deviation is seen at the low frequency 
end of the correlation. Ignoring the fiat-top noise, whose 
centroid frequency is zero, assuming Lp represents the 
fundamental QPO, the sequence of the peaked compo- 
nents is therefore 0.5-1-1.5-2. Notably, if we take the 
sub-harmonic peak Ls as the fundamental, the sequence 
becomes 1-2-3-4. In observations 6 through 11, all with 
vp around 3 Hz, an additional peak is detected, in three 
of these observations only at a 3cr level. Its frequency is 
consistent with being Zvp. For these observations, the 
harmonic sequence would then be 1-2-3-4-6. In Fig. [2l 
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Figure 1. Examples of the two types of PDS with type-C QPO. 
(a) PDS from one of the first 5 observations (Obs. 30188-06-01-01). 
Tlie solid line shows the best fit, while the dashed line represents 
the separate Lorentzian components (see text), (b) PDS of an ob- 
servation after the first five (Obs.: 30191-01-19-00). In both panels, 
the solid line shows the best fit, while the dashed line represents 
the separate Lorentzian components (see text). 



some components show systematic deviations from the 
harmonic lines. At low frequencies they correspond to 
Lpn and Lft, i.e. to the broad components in the fit. 
While L ft is not proposed here to be harmonically related 
to the others, the deviation of Lp„ must be discussed. 
The component is broad and difficult to characterize. In 
particular, the exactness of the Lorentzian approxima- 
tion is difficult to evaluate for such component (see Fig. 
[T]). This is more the case for the first observations, when 
the components is at lower frequencies. These deviations 
could be due to these effects, but their nature and pres- 
ence must be investigated in other sources. 

3.2. Frequency evolution 

Although a considerable amount of data showing type- 
C QPOs have been collected with RXTE, the nature 
of the frequency variations remains unclear. During 
the initial rise of the outburst of XTE J1550-564, on 
average, the type-C QPO between subsequent obser- 
vations rose monotonically with time and with count 
rate. [Sobczak ct al. (2000a) showed that the LFQPO fre- 
quency is directly correlated with the amount of (unab- 
sorbed) disk flux observed in the 2^20 keV band. There- 
fore, we first derived the disk count rate contribution for 
each observation, using the total count rate per PCU 
no rmalized by disk flux r atio using the best fit models 
by ISobczak et all (|2000b| l. assuming the ratio remains 
constant within one observation. We then produced a 
frequency vs. disk count rate relation (Fig. A linear 
relation is evident from the figure. For all observations, 
we produced a similar plot, but with a time resolution 
of 128 s (the disk rate ratio was assumed to be constant 
during the observation). The inset in Fig. |3] shows that 
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Figure 2. Correlations between characteristic frequencies of all 
Lorentzian components detected in the HIMS as a function of the 
frequency up of the Lp component. Circles, squares, filled circles, 
crosses and diamonds represent the second harmonic u^, peaked 
noise Vpn, fundamental vp, sub-harmonic Us, flat-top noise i/ft 
frequencies respectively (see Fig. [TJ. The dashed lines represent 
linear correlations with factors 2, 1.5, 1, 1/2 and 1/5 from top to 
bottom. 



the short-time resolution points (limited to one exam- 
ple observation) follow a similar correlation. The linear 
coefficients of the two correlations (long and short time 
scale) are inconsistent, most likely because of the choice 
of a constant disk rate ratio. Notice that correcting the 
single QPO peaks at 128 s resolution are not compatible 
with being a coherent oscillation. 

3.3. QPO width and coherence factor 

It is important to compare the widths of the different 
QPO peaks in order to uncover the nature of the quasi- 
periodic signal. In Fig. 21 we plot the FWHM of the 
Lorentzian components vs. their frequency for all the 
peaked components in the PDS. For comparison, for the 
zero-centered Lft component, we plot their FWHM vs. 
1/10 of j/p. It is clear that for all peaks, the FWHM is 
positively correlated with the frequency. The lines in Fig. 
0] indicate constant quality factors Q = fo/A. All peaks 
change their frequency while maintaining approximately 
the same Q. Lp and Lh have the same Q (around 11), 
while Lg is consistently broader (Q^3). The peaked noise 
component Lp„ is by definition even broader, but its Q is 
also consistent with a constant, 0.55, strengthening the 
identity of this component. Notice that for the peaked 
components, the characteristic frequency practically co- 
incides with the centroid frequency; the only deviation 
from this statement is of course the very broad Lp„. 

The component Ls (subharmonic) is broader than the 
other peaked components. This cannot be due to fre- 
quency modulation, as a different amount of modulation 
would destroy the harmonic relation on short time scales 
(the harmonic relation should hold on all time scales). 



4 Rao et al. 



7 




100 200 300 400 500 600 700 800 

Disk Count rate per PCU (s"') 

Figure 3. Fundamental QPO Lp frequency vs. disk count rate. 
Each cross corresponds to one observation. The inner panel shows 
the same relation on 128 s resolution from observation 30191-01- 
17-00. A typical error bar in frequency is shown in the inner panel. 



Frequency (Hz) 

Figure 4. Width vs. frequency relation for the peaked Lorentzian 
components in the PDS, plus the width of the Ljt component 
plotted as a function of 1/lOth of the Lp frequency. Symbols are 
the same as in Fig. [2] The dashed lines indicate Q factors of 0.55, 
3 and 11 from top to bottom respectively. 



It is however possible that this component undergoes an 
additional amplitude modulation and the L ft component 
is the prime candidate for providing it. In order to test 
it, we plot the correlation between the widths of the Ls 
and the Lft components in Fig. [S] A weak positive 
correlation is visible. The correlation coefficient is 0.68, 
corresponding to a chance probability of 1.4 x 10~^; a 
linear least-squares fit yields a value of 0.46. 

As the QPO fundamental frequency shows a good cor- 
relation with disk count rate on '^100 s timescales, it is 
possible that short-term variations of the centroid fre- 
quency due to this correlation are responsible for the 
width of the QPO peak. In order to estimate the amount 
of broadening due to this measurable frequency modula- 
tion, for each observation we first generated a 1-s binned 
light curve, then calculated its standard deviation. The 
resulting linear relation between the width oi Lp and this 
standard deviation is shown in Fig. [S) Notice however 
that the linear coefficient is much larger than that in Fig. 

m 

3.4. Energy dependence 

We performed an energy dependence analysis for all 
observations with a type-C QPO. Also type-B QPOs 
were analyzed for comparison. The two types of QPO 
(and harmonics) sho w different ene r gy de pendences, be- 
yond what shown bv iCasella et al.l ()2005() . An example 
is shown in Fig. [7] (this relation is similar for all other 
observations). In type-C observations, the rms of the 
fundamental and sub-harmonic peaks increase with en- 
ergy, while that of the second harmonic increases below 
a few keV, then decreases at higher-energies. In type-B 
observations, the rms of the fundamental and harmonic 
peaks increase with energy, while it is the rms of the 
sub-harmonic that shows a non-monotonic behavior: it 



1.5r 



0.5- 




0.5 



1.5 



2 2.5 

i/,widtli (Hz) 



3.5 



Figure 5. Width of the sub-harmonic peak Ls vs. width of the 
flat- top noise Lfi component. 



increases below 10 keV and decreases above. A similar 
relation was shown in Figs. 4 and 5 of ICui et al 



1999 ) 



for type-C QPOs and in Fig. 16_ofjHoman et aL| (pOOlf ) 
for type-B QPOs. Notice that iCui e t al. (199^ used a 
slightly different model for the power continuum. 

3.5. Interplay of components 
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Figure 6. Width of the fundamental QPO component Lp vs. 
standard deviation of disk count rate. Each point represents one 
observation. 
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Figure 8. Rms ratio (see text) vs. Lp frequency relation for both 
harmonic and sub-harmonic QPO. 
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Figure 7. Energy dependence for typical examples of type-C (left; 
Obs.: 30191-01-14-00) and type-B (right; Obs.: 40401-01-58-00). 
(a,b) PDS for high (13-33 keV; up) and low (2-3 keV; bottom) 
energies. The solid lines indicate the best fit; In (b), the high- 
energy PDS is shifted upwards by a factor of 10 for clarity. (c,d) 
rms dependence of the QPO peaks. The filled circles, open circles 
and crosses represent the Lp, Lf^ and Ls components respectively. 



3.5.1. QPOs 

In the past, the main focus has gone to the fundamental 
peak (Lf), while the harmonic peaks have received less 
attention (see e.g. iSobczak et al... 2000a for XTE J1550- 



564). We explored the interplay between components, 
i.e. how they vary differently both with energy and be- 
tween observations. The energy dependence described 
in the previous section indicates an interplay, where the 
harmonic rms ratio decreases with energy. This indicates 
that the signal becomes more sinusoidal at high energies. 
We define rms ratio as the rms of the harmonic (or sub- 
harmonic) peak divided by that of the fundamental. In 
Fig. IS] we plot the rms ratio vs. the Lp frequency. It is 
clear that the harmonic ratio decreases with the Lp fre- 
quency, while the sub-harmonic ratio increases, although 
with large scatter. 

3.5.2. Broadband components 

The same rms-ratio analysis can be performed with 
broadband components. The Lft rms ratio is defined as 
the rms of Lft divided by that of Lp„. We show their 
dependence on energy and Lp frequency in Fig. [S] The 
L ft ratio increases both with L p frequency and with en- 
ergy above 10 keV. Comparing the top panels in Fig. 
[9l one can see that the overall appearance of the PDS 
at low and high energies is quite different, as can also 
be seen in Fig. 4 of Cui et al. (1999). At low ener- 
gies, Lp„ is stronger and a break is visible around the 
QPO peak; at higher energies, Lft dominates and the 
break can be s een at low er frequencies, consistent with 
the Wiinands fc van der Klis (1999) correlation. This re- 
sult explains the presence of t he second (parallel) corre- 
lation in additio n to that of Wiinands fc van der Khal 
(|1999( ). found bv iBelloni et all (|2002| ). where the char- 
acteristic frequency of a broad-component is around the 
QPO frequency. 

4. DISCUSSION 

Following IBelloni et al.l (|2002D . we decomposed the 
PDS in Lorentzian components and found that two dif- 
ferent models are needed for the LHS observations (first 
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Figure 9. (a) best-fit models to power spectra of observation 
30191-01-14-00 (solid line: 13-33 keV; dashed line: 2-3 keV). (b) 
best-fit models to power spectra of observations 30191-01-23-00 
(solid line; QPO at 5.6 Hz) and 30188-06-05-00 (dashed line; QPO 
at 2.93 Hz), (c) energy dependence of the Ljt rms ratio (see text) 
for observation 30191-01-11-00. (d) Lj^ rms ratio versus Lp fre- 
quency for the full set of observations. 



five) and the HIMS ones. After investigating tlie inter- 
play among components, we argue tliat tlie transition is 
a smootfi process. As the disk flux increases, the Lp 
frequency increases, the Lbln ratio decreases and then 
disappear, the Lft and Ls ratio increase, the Lp„ and 
Lh ratio decrease. A similar interplay takes place as the 
energy increase. We can also conclude that Lp„ and Lh 
dominate at low Lp frequency, while Lft and Lg domi- 
nate at high Lp frequency. All these variations indicate 
that, depending on energy and time of the observation, 
the dominant low-frequency component might be Lpn 
or Lft, explaining the parallel correlation reported by 
IBelloni et all (12003) . 

Since we don't know the nature of the components, it's 
difficult to establish the physics of the interplay. How- 
ever, if we only consider the harmonic ratio, it is clear 
that the ratio is smaller at high energies and high disk 
flux, indicating that the QPO signal becomes more sinu- 
soidal. 

We studied the QPO components independently, with- 
out assumptions on their harmonic relation. We found 
that if we include the (broader) Lpn component, the fre- 
quencies of the four components Ls, Lp, Lpn, Lh follow 
a 1:2:3:4 sequence. This suggests that the Lg compo- 
nent, usually referred to as "sub-harmonic" , represents 
the fundamental frequency of the quasi-periodic oscilla- 
tion. In a few observations, an additional narrow peak at 
ivp was found. This corresponds to 6 times the subhar- 
monic frequencies, with no peak at 5 times is detected. 
Notice that there is no priori reason to expect that all 
harmonics are present in the signal at a detectable inten- 
sity. In other words, the absence of a higher harmonic 
does not weaken the 1-2-3-4 sequence, while the pres- 



ence of additional peaks at 0.5 and 1.5 times the funda- 
mental are difficult to reconcile with a harmonic series. 
Of course the detection of a peak at 2.'bvp would have 
strengthened the case. We also showed that the interplay 
between components can mask the presence of such a Lg 
peak. One problem with this interpretation is that Lpn 
is much less peaked than the other QPO components, 
with a Q factor of ~0.55 (see Fig. Moreover, at low 
frequencies, a deviation is observed from 3 times the fun- 
damental. As Fig. U shows, this additional broadening 
cannot be even larger than what could be caused by Lft- 
A component has been observed to change from being 
a clear band- limited noise compo nent into a clear QPO 
or vice versa ()Di Salvo et alll2001[ ). but it is still difficult 
to explain the difference between Lpn and other QPOs 
(see below). Alternatively, adopting Lp as fundamen- 
tal, two peaks would be present at unexpected frequen- 
cies: 0.5 (sub-harmonic) and 1.5 times the fundamental. 
Sub-harmonic signals are observed in nature, but they 
are associated to non-linear resonance, therefore involv- 
ing multiple frequencies, adding complication that linear 
models do not have. As a simple example, an oscilla- 
tor subject to a periodic force can develop sub-harmonic 
modes (at ujq/ti) if the period of t he force is a multipl e 
of that of the oscillation (see e.g. IButikovl[2002L [20081) . 
However, in order to overcome the the effects of damping, 
energy must be transferred to this mode through non- 
linear coupling. The presence of a 1.5i^i^ feature would 
imply another resonance, adding further complications. 
A simple harmonic series is a linear decomposition of 
a non-sinusoidal periodic component, without the need 
of involving non-linear effects. Since its first detections, 
the possibility that Vs represents the fundamental fre- 
quency, while sometim es not b eing detected, has been 
considered (see .T akizaw a et al.l [r997: Bclloni et al. 199'^; 
iRemillard et al.l l2002f ) . However, since no detection of 
a narrow peak at l.bvp was found, this assumption, al- 
though it could not be excluded, could not be proven. 

All frequencies are observed to increase with the count 
rate associated to the disk component. This is also true 
on short (128 seconds) time scales, although with dif- 
ferent parameters. A similar relation has been observed 
in typ e-C QPO of GRS 1915+105 also on sho rt time 
scales (jMarkwardt et al.lll999l: iMuno et al.lll999[) . GRS 
1915-1-105 is a very peculiar obje ct when observed at time 
scales longer than a second fsee lFender fc Bellonil 120(3^ 
and we can show now that this is also true for another 
less extreme black- hole transient. The difference between 
long- and short-time scale correlations could be intrinsic 
or just caused by error estimation of the disk flux ratio. 
If the latter is true, the correlation could in principle be 
used to guide the spectral fitting at short time scales. 

The result of the comparison between the FWHM of 
different QPO peaks is suggestive. The two simplest 
models for the broadening that gives origin to a QPO 
peak are amplitude modulation (AM) and frequency 
modulation (FM). Although these are indistinguishable 
in the PDS where only one peak is observed, when mul- 
tiple peaks are present, if the signals are subject to the 
same modulation, the effect of FM and AM is different. 
In AM, since the broadening is caused by the windowing 
due to the varying amplitude, the absolute width of the 
different peaks will be the same. In FM, the frequency 
modulation is multiplicative: this translates into peaks 
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that have the same Q factor (see e.g. Ivan der Klisl[T99l 
iBelloni et al.l 119971 ). Our results are complex. The Lp 

and Lh components have on average the same Q (no- 
tice that the scatter of individual measurements is large, 
but the overall correlation appears tight and consistent 
with an equal Q). This suggests that these two signals 
are frequency modulated. However, the "sub-harmonic" 
component Lg is considerably broader and the Lp„ at 1.5 
times the frequency of the fundamental is even broader. 
These two components must be subject to additional 
broadening. It is difficult (although it cannot be ex- 
cluded) to postulate an additional frequency broadening, 
since this would break the harmonic relation between fre- 
quencies on short time scales. It is more reasonable to 
assume that these two components are broadened fur- 
ther by an additional amplitude modulation, which af- 
fects only these two components, possibly associated to 
the Lft noise (see Fig. [2]). Not only we observe that 
Lp and Lh have the same Q factor. They maintain it 
throughout the outburst, which implies the modulation 
effects do not vary with time. 

We note that there is a QPO in a black-hole binary 
which has been analyzed in the time domain, due to 
strong statistics. It w as observed in GRS 1915-fl05 and 
it probably of type-C ([Morgan et al.lll997f) . The authors 
find that the oscillation shows a random walk in phase. 
Our result is consistent with this, as frequency and phase 
are directly related: frequency modulation is a type of 
phase modulation. Unfortunately, our signal is not suf- 
ficiently strong to allow such a detailed analysis in the 
time domain. 



5. CONCLUSION 

We have analyzed 47 RXTE observations with type C 
LFQPO acquired during the first half of XTE J1550-564 
1998 outburst. Satisfactory fits to the power spectra of 
the 6th-47th observations were obtained with a model 
consisting of flat-top noise Lft, peaked noise Lp„, fun- 
damental QPO Lf, sub-harmonic QPO L^, second har- 
monic QPO Lh and sometimes a third harmonic QPO 
components. We identify the peaked noise Lpn as a new 
harmonic component at l.bv the fundamental. This sug- 
gests that what we called previously sub-harmonic quasi- 
periodic oscillation (QPO) may actually be the funda- 
mental. 

A similar Q-factor between the fundamental Lp and 
the harmonic Lh was observed, suggesting a frequency 
modulation as cause of their width, whereas the sub- 
harmonic Lg is broader, and the I.Sz/q feature Lp„ is 
even broader. 

We also found a significant interplay among both 
QPOs and broad-band components with both frequency 
and energy. As the disk flux increases, Lp frequency 
increases, Lft and Ls ratio increase, while Lp„ and Lh 
ratio decrease. A similar interplay also takes place as 
the energy increase. The nature of the interplay remains 
unknown, since we didn't even know the nature of the 
components. However, from the interplay we may ex- 
plain the PDS difference in the first 5 observations, and 
suggest the transition from LS to HIMS is a smooth pro- 
cess. 
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